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We present a summary of our results for the light hadron spectrum in quenched lattice QCD 
and preliminary ones in two-flavor full QCD. For the quenched hadron spectrum we find that the 
mass formulae predicted by quenched chiral perturbation theory(QChPT) give a good description 
of our results. Employing the QChPT mass formulae for chiral extrapolations we conclude that the 
mass spectrum in the continuum limit deviates unambiguously and systematically from experiment. 
For our two-flavor full QCD results we focus on the dynamical quark effects on the light hadron 
spectrum. 



I. INTRODUCTION 

Numerical studies of QCD on the lattice provide us a quantitative understanding of the dynamics of strong interac- 
tions. One of the main goals is to derive the hadron mass spectrum from the first principle. In order to give a precise 
prediction with lattice QCD simulations we have to control the following three physical sources of systematic errors 
besides the quenching error: (i) scaling violation, (ii) finite size effects and (iii) chiral extrapolations. 

In the quenched approximation the GFll collaboration carried out an extensive calculation to reduce these three 
systematic errors and found that the quenched spectrum is consistent with experiment within 5 — 10% errors Q. The 
aim of our work is to tell definitely how the quenched spectrum deviates from experiment, by considerably diminishing 
the statistical and systematic uncertainties. 

On the other hand, in spite of much effort devoted in the full QCD spectrum study we have not yet been able 
to answer a fundamental question: to what extent the dynamical quarks affect the light hadron spectrum. Before 
embarking upon a realistic full QCD calculation in the future we attempt to settle this question. 

In the first part of this report we present our results for the quenched light hadron spectrum, examining the validity 
of QChPT for the chiral extrapolation. The two-flavor full QCD results obtained so far are presented in the second 
part, where we investigate possible signs of dynamical quark effects on the light hadron spectrum. More details are 
found in Ref g. 

II. QUENCHED LIGHT HADRON SPECTRUM WITH THE WILSON QUARK ACTION 

A. Details of numerical simulation 



TABLE I. Simulation parameters for quenched QCD. 
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Our simulations are carried out with the plaquette gauge action and the Wilson quark action. To control the 
systematic errors we carefully choose our run parameters, which are summarized in Table ||. We employ four (3 values 
so that the lattice spacing covers the range a ~ 0.1 — 0.05fm to remove the scaling violation effects by extrapolation 
of the data to the continuum limit. To avoid finite size effects we keep the physical spatial lattice size approximately 
constant at La « 3fm. A previous study showed that the finite size effects are 2% or less already at La sa 2fm ||^ . The 
most subtle issue in controlling the systematic errors is associated with chiral extrapolations of the hadron masses. 
Although the mass formulae predicted by QChPT are considered to be plausible candidates for the fitting functions 
of chiral extrapolations, their validities should be checked employing a wide range of quark masses. At each value 
of j3 we choose five values of the hopping parameter K corresponding to m-^jmp « 0.75, 0.7, 0.6, 0.5 and 0.4. The 
heaviest two values m-^lrfip « 0.75 and 0.7 are taken to be around the physical strange quark mass. In terms of these 
quark masses we calculate hadron masses both for degenerate and non-degenerate cases. 

Gauge configurations are generated with the 5-hit pseudo heat-bath algorithm incorporating the over-relaxation 
procedure quadruply and employing periodic boundary condition. Quark propagators are solved in the Coulomb 
gauge both for the point and the smeared sources with periodic boundary conditions imposed in all four directions. 
For the hadron mass measurement we use the smeared source quark propagators. 

The physical point for the degenerate up and down quark mass is determined by to^ (135.0) and (768.4). For 
the strange quark mass we employ (497.7) or m^(1019.4). The lattice scale a~^ is set with Wp. 



B. Chiral extrapolations 



We first examine the validity of QChPT investigating the presence of quenched chiral logarithm in the pseu- 
doscalar(PS) meson sector. For the quark mass dependence of the PS meson mass mpg the QChPT predicts W 
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where m and rUs are masses of two valence quarks in the PS meson. This ratio diverges logarithmically toward the 
chiral limit due to the S term. To detect the contribution of the S term in a direct manner we introduce the two 
variables: x ^ 2 — (nis + m)ln(ms/rn)/(ms — rn) and y ~ Amms/ {mg + m)^ x / {m^m^^) , for which cq.(|l|) leads 
to the relation y = \ + 5 ■ x. Here tt and r/ are the degenerate PS mesons with quark mass m or and K is the 
non-degenerate one with m and m^. It should be noted that the quark masses are defined by an extended axial 
vector current Ward identity (AWI) V^A)f* = 2mf^^P |], where we are free from ambiguities originating from the 
determination of the critical hopping parameter Kc- Figure |^ shows the distribution of y as a function of x. Our data 
fall in a wedge shaped with y = 1 + 0.08x and t/ = 1 + 0.12.T. A PS meson decay constant ratio y = ff(/{fnfij), for 
which QChPT predicts y ~ 1 — 6/2 ■ x, is another quantitative test for 6. Our data indicate that the value of 6 is 
within 6 = 0.08 - 0.16. 

Since we observe clear evidence for the existence of the quenched chiral logarithm, we now try to fit the PS meson 
masses with the functional form of eq.(|l|). In Fig. ||we show a typical result for chiral extrapolations of the degenerate 
PS meson mass and the AWI quark mass. There are two features to be remarked. One is a good description of the 
AWI quark mass with a linear function of l/K, which means the AWI quark mass is proportional to the vector Ward 
identity(VWI) quark mass m^^^ = {l/K — l/Kc)/2. The other is a good agreement between the critical hopping 
parameter Kc determined from the QChPT fit of mp^ and that from the linear fit of m^^^. They are consistent 
within 2.5cr, which should be compared with 17fT (12fT) discrepancy for the case of a quadratic (cubic) polynomial fit 
of mpg. 

Our examinations for the chiral properties of the PS meson masses strongly suggest the validity of QChPT with 
6 « 0.10(2). It would be legitimate to apply the QChPT mass formulae for the chiral extrapolation of the vector 
meson masses and the baryon ones. 
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FIG. 1. Test of the quenched chiral logarithm using 
the PS meson masses. 
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FIG. 2. Chiral extrapolations of the degenerate PS 
meson mass and the AWI quark mass as a function of 
XjK 5.9. 



For the vector meson masses we make simultaneous but uncorrelated fits of the degenerate and non-degenerate data 
employing the functional form predicted by QChPT: 
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where C1/2 — —^T^g^^ with a phenomenological coupling constant of the vector meson quenched chiral Lagrangian. 
It should be noticed that the Oiraps) term is a characteristic of QChPT. This model function up to the 0{mpg) 
terms describes our vector meson mass data well as illustrated in Fig. ||. A small bending toward the chiral limit 
found in the fitting result of mp reflects the contribution of the 0{mps) term. However, the fitted values for the 
coefl[icient C1/2 of 0{mps) term are much smaller than expected. The average value of C1/2 over four /3 points gives 
Ci/2 = —0.071(8), whose magnitude is ten times smaller than the phenomenological estimate C1/2 — — 47rf/|5 w —0.71 
with 5 = 0.1 and 32 — 0.75. The magnitude of the 0{mps) term in eq.(^) is 0.07 x TOtt ~ lOMeV, which means about 
a 1% contribution to mp. If we employ the fitting function up to the 0{nVpg) terms, we find a few times larger value 
for Ci/2 at each beta. In this case, however, the fitting results are very unstable against the value. 
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FIG. 3. Chiral extrapolations of vector meson masses 
as a function of at /3 = 5.9. Open symbols in the inset 
represent extrapolated values at the physical degenerate 
up and down quark mass. 
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FIG. 4. Same as Fig. ^for octet and decuplet baryon 
masses. 
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FIG. 5. Continuum extrapolations of quenched 
hadron masses with rriK as input for the strange quark 
mass. 
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FIG. 6. Quenched light hadron spectrum in the con- 
tinuum limit. GFll results [Q] are also plotted for com- 
parison. Horizontal bars denote experimental values. 
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For fitting of the octet and decuplet baryon masses we take the same strategy as for the vector meson mass case. 
Assuming the QChPT mass formulae we fit the degenerate and non-degenerate data simultaneously without including 
correlations. While we find that 0{mps) and 0{'m^g) terms are sufficient to reproduce our decuplet baryon mass 
data as shown in Fig. ^, for the octet baryon masses a O(mpg) term is inevitably required to explain the convex 
curvature of the nuclcon mass data. The fitted value for the coefficient C1/2 of the 0{mps) term in the nucleon is 
Ci/2 = —0.118(4) as an averaged value over four f3 points. This is less than half of the phenomenological estimate 
Ci/2 = -(37r/2)(i:» - 3F)2<5 « -0.27 with 6 = 0.1, F = 0.5 and D = 0.75. The Oimps) term contribution to the 
nucleon mass is about 16MeV or 2%. We find a similar value C1/2 — —0.14(1) for A. If the 0{m'^g) terms are 
included in the fitting function of the decuplet baryon masses, we obtain much smaller value for C1/2 in A. 

Whereas we observe that the mass formulae based on the QChPT reproduce our data adequately, it is instructive 
to try another chiral ansatz employing polynomial fitting functions in 1/K (cubic for N, quadratic for others). We 
employ independent fits for the different degeneracies in the vector meson and the octet and decuplet baryons. The 
fitting results are shown by dashed lines in Figs. H, ^ and|4[ Although the QChPT fits seem to give better descriptions 
of our hadron mass data than the polynomial ones, the differences between them are fairly small. 



C. Quenched light hadron spectrum in the continuum limit 

The lattice spacing dependences of the hadron masses are shown in Fig. ^. We expect that the leading scaling 
violation effects are of 0{a) in the Wilson quark case. A linear extrapolation to the continuum limit assuming the 
fitting function m — mo{l + aa) yields a ^ 0.2GeV, for which we find that higher order terms are safely negligible, 
e.g., (aa)2 - 0.01 at a = 0.5GeV-i. 

While the results of the QChPT fits and the polynomial ones differ by about 3%(5cr) in the largest case after 
the chiral extrapolation at each (3, the differences in the continuum limit are within 1.5%(1.5(t) of the results of the 
QChPT fits. At a few percent level of statistical errors it is hard to appreciate the differences between the QChPT 
chiral extrapolations and the polynomial ones. 

Our final result of the quenched light hadron spectrum is presented in Fig. ^ The spectrum deviates from experi- 
ment systematically and unambiguously. The following discrepancies should be noticed: The K-K* meson hyperfine 
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splitting is definitely underestimated by 9.5%(4.3(7) for the ttik input case and by 16%(6.1(t) for the input case; 
The nucleon mass is appreciably smaller than experiment by 7%(2.5a); For the strange octet baryon masses our data 
with the rriK input are systematically smaller than experiment by 6 — 9% (4 — 7 a), which is much reduced in the 
input case; While the tjik input leads 30% smaller estimates for the decuplet mass splittings in average, we observe 
rather good agreement for the input results. 



III. LIGHT HADRON SPECTRUM IN TWO-FLAVOR FULL QCD 
A. Details of numerical simulation 

The simulation of full QCD requires a huge amount of computing time compared to the quenched approximation. In 
order to secure an adequate physical spatial lattice size for avoiding finite size effects, the lattice spacing is compelled 
to be coarse a~^<2GeV, which urges us to employ improved actions. Our choice is a RG-improved gauge action ||] 
and the SW quark action |Q with a mean- field improved value of csw- This decision is based on a former comparative 
study of various combinations of improved gauge and quark actions at a^^ « IGeV ||^. 

Our full QCD study is performed with two flavors of sea quarks which are supposed to be the degenerate dynamical 
up and down quarks. Our simulation parameters are summarized in Table 0. Four (3 values covering the lattice 
spacing in the range a^^ w 1 — 2GeV are employed to examine the scaling behavior of the light hadron spectrum. 
We keep the physical spatial lattice size approximately constant at La 2.4fm except for the finest lattice where 
La « 2.0fm. At each value of f3 we choose four values of sea quark mass corresponding to niT^/mp w 0.8, 0.75, 0.7 and 
0.6. Gauge configurations including the dynamical sea quark effects are generated with the HMC algorithm. 

We calculate hadron masses using the five valence quarks whose masses correspond to m-^/mp « 0.8, 0.75, 0.7, 0.6 
and 0.5. Hadron propagators are constructed with the smeared source quark propagators at every fifth trajectory. 
Hadron masses are extracted by an exponential fit ignoring correlations between time slices. Errors are estimated by 
the jackknife method with a bin size of 10 configurations (50 HMC trajectories). 

The light hadron spectrum is obtained by setting the sea quark mass on the physical up and down quark mass 
and the valence quark mass on the physical up and down quark mass or the strange one. The physical point of the 
degenerate up and down quark mass is fixed with TOtt and nip as input. For the strange quark mass the physical point 
is determined from nix or m^. The lattice scale is set with rup. 



B. Chiral extrapolations in the sea and valence quark masses 

In Fig. 1^ we plot a result for the sea and valence quark mass dependences of hadron masses at j3 — 1.95. The 
data are parameterized by the averaged hopping parameter 1/Kav = (l/^t)a((i) + l/^'i)a((2))/2 of the two kind of 
valence quarks constituting the hadrons. S represents a valence quark with Kyai = K^ea and V a valence quark with 



TABLE II. Simulation parameters for full QCD. 
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We observe that the data of the meson masses are distributed on different four hues in accordance with the four 
sea quarks. Partially quenched data on each sea quark show almost linear behavior both for m^g and mvec- Their 
slopes, however, slightly depend on the sea quark masses: As the sea quark mass decreases the slope decreases for 
TOpg, while the slope increases for myec- In the S-S channel nr^g behaves almost linearly in 1/Ksea, whereas mvec 
exhibit a convex curvature. Chiral extrapolations of the meson masses are made by global fits of all the data assuming 
quadratic functions in terms of l/Kgea and l/Kyai- For the PS meson mass we employ 
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where bare quark masses are defined by mg^g^fya.i{i) = (l/'^sea/t)ai(i) ~ l/-f^c)/2 and rn^ai is the averaged value of two 
valence quark masses. Similar fitting function without the valence- valence cross term 'Tii,^/(]^)m„(j;(2) is used for the 
vector meson masses and the decuplet baryon masses. The fitting results for the meson masses are drawn in Fig. 

As for the octet baryon masses we find a rather complicated situation in Fig. Partially quenched baryon masses 
are not functions of m^o;. For the fitting function we take a combination of linear terms based on the ChPT prediction 
and general quadratic terms of individual myai(i)- This function has 12 free parameters in all, which are determined 
by a combined fit of E-likc and A-likc baryon masses. 
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FIG. 7. Sea and valence quark mass dependence of 
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C. Sea quark effects on the light hadron spectrum 



The results for the light hadron spectrum are presented in Fig. 1^ as a function of the lattice spacing. For comparative 
purpose we also plot the quenched results(open symbols) from the last section, which would help us distinguish the 
dynamical sea quark effects. 

We observe an intriguing feature in the meson sector: The K* and (p meson masses in full QCD are heavier than 
those in quenched QCD at finite lattice spacing and become closer to experimental values as the lattice spacing 
decreases. However, our full QCD results do not necessarily agree with experiment in the continuum limit. We have 
two possible reasons. One is the relatively heavy sea quark masses corresponding to TOjr/^p ~ 0.6 — 0.8, for which 
the chiral extrapolations are rather ambiguous. The other is the quenched treatment of the strange quark. 

For the baryon sector it is hard to read any meaningful implication in comparison between the full QCD results and 
the quenched ones. At the finite lattice spacing the full QCD spectrum is consistent with the quenched ones within 
rather large errors, while the lattice spacing dependence of the decuplet baryon masses in full QCD is opposite to 
that in quenched case. 



IV. CONCLUSIONS 

We have presented our results for the quenched light hadron spectrum. After examining the validity of QChPT 
in the pseudoscalar meson masses we applied its hadron mass formulae to the chiral extrapolations. Although the 
QChPT predictions reproduce our hadron mass data well, for the vector meson and baryon masses it is hard to confirm 
the contributions of the term specific to QChPT. The spectrum in the continuum limit deviates systematically and 
unambiguously from experiment on a level of about 10% far beyond the statistical error of 1 — 2% for mesons and 
2 — 3% for baryons. 

Our full QCD study is an exploring step toward a realistic QCD simulation. For the light hadron spectrum we 
found an encouraging result that the deficiency of the meson hyperfine splitting in the quenched approximation is 
largely compensated with the dynamical sea quark effects. At this stage, however, the most important conclusion 
is as follows: To perform a close investigation of the sea quark effects we need a direct comparison between the full 
QCD results and the quenched ones employing the same gauge and quark actions with the same a^^, mT^/nip and 
La. Work in this direction is now in progress. 
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